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Flat galaxy rotation curves and the accelerating Universe both imply the existence of a critical accelera- 
tion, which is of the same order of magnitude in both the cases, in spite of the galactic and cosmic length 
scales being vastly different. Yet, it is customary to explain galactic acceleration by invoking gravitation- 
ally bound dark matter, and cosmic acceleration by invoking a 'repulsive' dark energy. Instead, might it 
j— 1 not be the case that the flatness of rotation curves and the acceleration of the Universe have a common 

cause? In this essay we propose a modified theory of gravity. By applying the theory on galactic scales 
we demonstrate flat rotation curves without dark matter, and by applying it on cosmological scales we 
demonstrate cosmic acceleration without dark energy. 
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Dark matter has been postulated to explain the flat non-Keplerian behaviour of rotation curves of galaxies, 
and the high dispersion velocities of galaxy clusters. Dark energy/cosmological constant has been postu- 
lated as an explanation for the observed late-time cosmic acceleration. Together, cold dark matter and a 
cosmological constant provide the standard AC DM model of structure formation and large-scale-structure 
which best fits the current cosmological data. Nonetheless, there are chinks in the armour! Apart from 
the well-known facts that dark matter has not yet been detected in the laboratory, and that the inferred 
value of the cosmolgical constant is too low compared to the theoretically preferred value by some 10 120 
orders of magnitude, there is another remarkable observation which the ACDM model does not seem to 
account for. And that is the existence of a universal acceleration cHq ~ 1CP 10 cm sec -2 , where Hq is the 
present value of the Hubble parameter, which prevails in the outer regions of galaxies, in galaxy clusters, 
and at the edge of the observed Universe: 
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In this essay, we seek an understanding of the cause for this universal acceleration, and for flat galaxy 
rotation curves and cosmic acceleration, by postulating a fourth order modified gravity as a common 
explanation, as an alternative to dark matter and dark energy. We propose, for reasons to be discussed 
below, that the current structure and evolution of the Universe is described by the following effective field 
equations 

BT - -g» u R = ^f-T^ + L 2 R^., (2) 
2 c ' 



The length parameter L is scale-dependent and defined as being proportional to the square-root of the 
mass M of the system under study, in such a way that GM/L 2 is a universal constant equal to cHq. This 
holographic assumption [mass proportional to area rather than volume] seems to be universally valid for 
large astrophysical systems (see e.g. Fig. 7 in [1]) and entirely consistent with the observation in Eqn. 
Q. 

We will solve these modified gravity equations for two different spacetime metrics: 

(i) The Newtonian weak-field non-relativistic approximation for the scalar potential 4>: 

„2 _ / -, , 2 <M 2 ,.2 ,2 i .2 > 2 



ds z = ( 1 + J c z dt z - dx z - dy z - dz l ; (3) 



and 



(ii) The spatially flat Robertson- Walker cosmological metric 

ds 2 = c 2 dt 2 -a 2 (t) [dx 2 + dy 2 + dz 2 } . (4) 

In the first case we will show how the resulting fourth order biharmonic modification of the Poisson equation 
explains flat galaxy rotation curves without dark matter, and in the second case we will show how a fourth 
order modification of the Friedmann equations explains late-time cosmic acceleration without dark energy. 

It can be shown that for the weak-field metric Q the modified equations ^ reduce to the following 
fourth order biharmonic correction to the Poisson equation [21 [3] 

V 2 <f> - L 2 V 4 <p = 4nGn(r) (5) 

where /z(r) is the matter density distribution of interest. From observations of galaxies this is known to be 

where 
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and 



MM = 4vr / dr'r u u(r') = M (7) 

Jo \r + r c J 

„ _ J 1 for HSB galaxies, , , 

P ~ \ 2 for LSB & Dwarf galaxies. ^ 

For this mass distribution, the biharmonic equation ^ can be solved to obtain the following remarkable 
Yukawa form for the radial component of the gravitational acceleration a = — V</> inside the galaxy 



GM(r) ( , M 

fl(r) = — 1 + V¥ 



1 - cxp( -/•/!)( l + j 
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where Mq is a constant of integration and M is the total mass of the galaxy [3]. The rotation curve v(r) 
of the galaxy is obtained by setting a(r) = v 2 {r)/r. If we choose L according to the prescribed relation 
GM/L 2 = cHq [this gives L ~ 10 kpc] and if we choose Mq = 10 12 Mq it is possible to fit the flat rotation 
curves of a large class of galaxies, without invoking dark matter [31 0]. It is to be noted that the values of 
L and Mq chosen here are consistent with all known laboratory and astronomical tests of departures from 
the inverse square law [5]. For r <C L, the inverse square law holds. Also, for r ^> L the inverse square 



law holds but with an effectively larger value of the gravitational constant. Interesting new physics arises 
around r > L. 

Greater insight into the modified acceleration obtained by us can be had by expanding the acceleration 
Q around r = L and writing it as a sum of two terms : a part that falls as 1/r 2 and is independent of L, 
and a part that falls as 1 jr and depends on L. This gives 




r(kpc) 

Figure 1: Red curve: velocity due to first term, Green curve: velocity due to second term, Blue curve: 
total velocity. 

The second term dominates for r > L. Fig. 1 above plots the rotation velocity curve, and it is clear 
that the Keplerian fall-off due to conventional Newtonian acceleration [first term] is modified to a flat 
rotation curve [second term] for r > L, because of the proposed generalisation of Einstein gravity. 
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Next, we turn to cosmology and write the modified Friedmann equations resulting from Eqn. ([2]), for 
the Robertson- Walker metric. Since one is no longer in the weak-field limit, the equations are complicated, 
but it is adequate for our purpose to write them in the form 



a 2 t2 _ , ...s 8ttG . . 

-^-L u J 7 i{a,a,a,a) = ——p (11) 



— + Tg- — L\jT%{a, a, d, a, "a) = 8-TrGp (12) 

where T\ and Ti are known functions of their arguments which we do not write explicitly. In accordance 
with the scaling principle, since we are now considering the entire observable Universe, we have Ljj = cHq 1 . 
It turns out that these equations possess a very interesting analogy with the weak-field case considered 
above. There is a space f-> time symmetry: the behaviour seen in the previous case around r > L is 
repeated here for t > Hq 1 . 

For times t <C Ljj the modifying gravity terms can be neglected and these equations reduce to the 
standard Friedmann equations which give a decelerating solution. This is analogous to the behaviour seen 
in the solution ^ for the gravitational acceleration: just as the inverse-square law holds for r <C L in 
Eqn. ([9|, for t -C Hq 1 one recovers the standard Friedmann evolution. Things start to get interesting for 
t > Ljj = Hq 1 . For the sake of simplification one can divide the evolution into two phases: (i) evolution 



for t <C Hq 1 , and (ii) evolution for t > Hq 1 . In the second case, one can show from Eqn. (12) for the dust 
case p = 0, assuming a power-law solution, that the evolution is dominated by the modifying terms, and 
the scale-factor is given by 

*> - (i3 » 

where T is the current age of the universe. One has an accelerating solution with a > 0. This is in complete 
analogy with the galactic case: the modifying term increases the acceleration, and in the present cosmo- 
logical case the rise is enough to change the sign from deceleration to postive acceleration. Furthermore, 
the acceleration decreases with increasing time, and it can be shown that for t 3> ^o" 1 the decelerating 
phase reappears. It is easily shown that T = 5-ff ( ^ 1 /2, and as can be expected the acceleration is of the 
order of cHq. [The vanishing of the scale factor at t = Ljj is only an approximation, and simply reflects 
the smallness of its value at t = Ljj compared to the present value a = 1; the true singularity of course 
occurs at the beginning of the previous decelerating phase] . 
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Using the solution (13) for the scale factor in the first Friedmann equation (11) it can be shown that 



the dust density evolves as 



p(t) oc i oc -ig. (14) 



The density falls with the scale factor slowly compared to the standard case [where p oc 1 /a 3 ] , because the 
modifying term in the Einstein equation implies that 



t + L 2 R^, ag ) = 0. (15) 



The onset of fourth order gravity at the present cosmological epoch, which we have proposed as a phe- 
nomenological representation of an underlying physics, implies an exchange of energy-momentum between 
gravity and matter fields at galactic and cosmological scales. This is one possible unified way of under- 
standing the otherwise unusual features of flat rotation curves and cosmic acceleration. 

The underlying physics may possibly have to do with the effect that averaging over small-scale inhomo- 
geneities has on dynamics on larger scales, a phenomenon which is perhaps not fully understood at present. 
In fact, such considerations led to the suggestion from Szekeres [6] of essentially the same equation as 
Eqn. ([2]) proposed here. 

Understanding structure formation in this picture is a challenge. But we have a unified description 
of critical acceleration on different scales, which comes, not from invoking different matter/energy com- 
ponents, but a phenomenological modification of gravity whose root cause may be the same at different 
scales, and which calls for further study. 
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